The compact core-jet region of the superluminal quasar 3C 216 
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Abstract 

Space very long baseUne interferometry (VLBI) observations of the quasar 3C 216 with the VLBI Space 
Observatory Programme (VSOP) reveal that the parsec scale structure of the source can be well described 
by compact jet models developed for interpreting the core region of radio- loud active galactic nuclei (AGN). 
The measured brightness temperature of Tb = 7.9 x 10^^ K is comparable to the inverse Compton limit, 
from which we determine a lower limit of S ^ 3.17 for the Doppler-boosting factor. The apparent transverse 
velocity of the superluminal component is /3app = (3.0 ± 0.2) assuming a constant velocity, but decel- 
eration of the jet material cannot be excluded from our data. A combination of the above values indicates 
that the viewing angle of the core-jet to the line-of-sight is less than 18?4, and the jet Lorentz-factor exceeds 
3.16. The observed small size of the source is probably caused by both interaction with the interstellar 
medium, and a projection effect. 
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1. Introduction 

The central region of radio-loud active galactic nuclei 
. (AGN) can be studied with sub-milliarcsecond (mas) res- 
olution using the space very long baseline interferometry 
(space VLBI) technique; this corresponds to only a few 
parsecs projected linear size even for the highest redshift 
sources. The inclusion of an orbiting element into the 
VLBI array is important because the achievable resolu- 
tion increases without using higher observing frequen- 
cies. As theory predicts that the size of the optically 
thick region of compact jets is inversely proportional to 
the frequency (Blandford & Konigl 1979), sources that 
are too compact for ground VLBI arrays can only be re- 
solved by space VLBI. The other advantage of this tech- 
nique over ground-based VLBI is that higher brightness 
temperatures can be measured. The maximum achiev- 
able brightness temperature (for an isotropic source) is 
limited by inverse Compton scattering to Tb = 10^^ K 
(Kellermann & Pauliny-Toth 1969); Tb values exceeding 
~ 10^^ K indicate Doppler-boosting of the emission. Ob- 



servational evidence for Doppler-boosting and investiga- 
tion of parsec-scale jet misalignments at sub-mas resolu- 
tion (e.g. Tingay et al. 1998) help us to refine orientation- 
based unification models. 

Below we present results of our VLBI Space Ob- 
servatory Programe (VSOP) observations of 3C 216 
(0906-f430). The source (z = 0.67, Spinrad et al. 1985) 
shows high optical polarization and variability, typical of 
the class of blazars (Angel & Stockman 1980). It was 
classified as a compact steep spectrum (CSS) quasar as 
well (e.g. Fanti et al. 1985; Pearson & Readhead 1988), 
though this classification may not be entirely correct as 
the observed small source size (~15 kpc) is probably not 
intrinsic, and within 100 mas of the centre the radio 
structure is reminiscent of core dominated sources (Fejes 
et al. 1992). There are slightly misaligned radio lobes 
on arcsecond scales oriented to South- West and North- 
East directions from the centre (Pearson et al. 1985; Fejes 
et al. 1992); the latter seems to align with the barely- 
resolved optical structure observed by the Hubble Space 
Telescope (de Vries at al. 1997). Images made by the 
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VLA at 1.4 GHz (Barthel et al. 1988) and by the MER- 
LIN array at 408 MHz (Fejes et al. 1992) reveal the pres- 
ence of a low surface brightness halo about 8" in size 
- there is no evidence for larger scale emission at lower 
frequencies, according to 81.5 MHz long baseline inter- 
ferometric observations (Hartas et al. 1983). Regard- 
less of the strong lobes that dominate the spectrum at 
low frequencies, 3C216 shows the structural properties 
of a core dominated source at centimeter wavelengths on 
VLBI scales. It has a bright core component and a jet 
that is perpendicular to to the arcsecond-scale radio lobes 
(Fejes et al. 1992). The jet follows a gently curved path, 
with underlying wiggles, and ends in a sharp bend at 
140 mas from the centre (Fejes et al. 1992; Akujor et al. 
1996). In the innermost part of the source, components 
emerge at superluminal speeds (Barthel et al. 1988). Ven- 
turi et al. (1993) showed that deceleration takes place in 
the inner few mas (10-20 parsecs). There is also a quasi- 
stationary component at ~ 1.5 mas from the radio core 
at 5 GHz. Our main goal was to investigate this region 
in more detail. 

Table 1. Characteristics of radio antennas used in our 
VSOP experiment. 



Name Diameter (m) SEFD^ (Jy) 



Effelsberg 100 20 

Westerbork 93*' 70 

Jodrell Bank 25 320 

Noto 32 260 

Onsala'^ 25 600 

Torun 32 250 

Medicina 32 296 

Green Bank 43 120 

HALCA 8 15300 



^ System Equivalent Flux Density. 

^ Used in phased array mode, an equivalent diameter is 
given. 

Did not observe due to high winds. 

2. Observations, data reduction and imaging 

The source was observed at 5 GHz with the array of 
the VSOP satellite HALCA, the 43m Green Bank tele- 
scope and the western part of the European VLBI Net- 
work (EVN) on 14/15 February 1999. Characteristics 
of the radio antennas are listed in Table 1. There were 
three tracking passes during the 8.5 hour experiment us- 
ing the Green Bank, Tidbinbilla and Goldstone tracking 
stations. The data were recorded in left circular polariza- 
tion in VLBA mode, and correlated at the NRAO VLBA 
correlator in Socorro (USA). We used AIPS (Cotton 
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Fig. 1.. Contour representation of the bending jet of 3C216, 
imaged using natural weighting of the data taken at 5 GHz. 
Due to the large image pixel size of one mas, only data from 
the shorter EVN baselines were included. Contour levels are 
-0.35, 0.35, 0.5, 1, 2, 4, 8, 16, 32, 64% of the peak bright- 
ness of 507 mjy/beam, the restoring beam is 8.3 X 6.3 mas at 
PA = 27°. The inset shows the wiggling jet reproduced with 
0.4 mas image cellsize and some uti-tapering, which resulted in 
a restoring beam of 3.0 X 1.4 mas at PA = —15° (about the 
resolution of the ground-only array of this experiment). Con- 
tour levels are -0.3, 0.3, 0.5, 1, 2, 5, 10, 25, 50, 99% of the 
peak brightness of 461 mJy/beam. 

1995; Diamond 1995) version 150CT98 for data process- 
ing at the FOMI Satellite Geodetic Observatory, Hun- 
gary. The measured antenna system temperatures (Tgys) 
were used for a-priori amplitude calibration where avail- 
able. We used nominal T^ys values for Green Bank and 
HALCA. Fringes were found with high signal-to-noise ra- 
tio for space baselines during the experiment, except in 
the short third tracking pass. We averaged data in fre- 
quency and time (one minute integration time) and per- 
formed imaging in DIFMAP (Shepherd et al. 1994). 



We experimented with a variety of uv cut-offs and data 
weighting in imaging, in order to image the source struc- 
ture on different spatial scales. We were able to recon- 
struct the jet bend at ~ 140 mas from the core using 
natural weighting - this feature was resolved out using 
uniform weighting, when ground-space baselines domi- 
nate in imaging. We found that the gain of the EVN 
telescopes - which have baselines sensitive to the 100 
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mas scale structure - could be better self-calibrated if 
we added clean components from the large scale emis- 
sion "manually", when performing imaging with uniform 
weighting. 

The process is described below. 1) A self-calibration 
and imaging cycle was carried out using natural weight- 
ing. 2) The source model file obtained was edited so that 
clean components within 5 mas from the phase centre 
were deleted (emission from larger scales do not have a 
significant contribution to the correlated flux density on 
the Earth-space baselines). 3) A new self-calibration and 
imaging cycle was started with uniform weighting of the 
data. After the starting phase self-calibration to a 1 Jy 
point source model, the edited model file from 1) was 
appended to the clean component model. Imaging of the 
inner few mas region of the source was then continued in 
the usual way. This method should be applied with cau- 
tion but proved to be useful in our case. While the basic 
structure of the source obtained as described above was 
not different from what can be achieved by the standard 
way, we were able to improve the dynamic range of the 
final space VLBI map. 

3. Results 

Contour representations of the resulting images - ac- 
cording to the various imaging strategies - are shown in 
Fig. 1. and Fig. 2. The bright knot at the jet bend can 
be identified at 140 mas from the core in Fig. 1., where 
the inset shows the wiggling jet in more detail. Compo- 
nents in the full resolution space VLBI image (Fig. 2a) 



are labelled in the same way as in Venturi et al. (1993). 
The core region is resolved; however, the clean beam is 
highly elongated roughly along the jet direction. The re- 
sults of model fitting of the self-calibrated data are listed 
in Table 2. Using a 0.5 mas circular restoring beam, 
Component A is apparently well separated from the core 
(Fig. 2b). The separation is estimated as 1.25±0.05 mas. 
This is not in agreement with the results of model fitting; 
an elongated component closer to the core (0.87 mas) 
gives the best fit to the data (see Table 2.). Component 
B is clearly resolved, and seems to be edge brightened 
on the western side. It is difficult to estimate the posi- 
tional uncertainty of component B; we will use 0.2 mas 
in further analysis (as well as for other epochs). The po- 
sition angle of the inner jet components {PA = 152°) is 
slightly misaligned with respect to the 100 mas scale jet 
{PA = 146- 149°). 

In order to scale angular sizes to projected linear sizes, 
we used Hq = 100 /ikms""'^ Mpc~^ and qq = 0.5; 1 mas 
corresponds to 3.9/i~^pc at the distance of 3C216. We 
also calculated the observed brightness temperatures of 
the source components in the core region, after correction 
for the z cosmological redshift, using the formula given 
by Kellermann and Owen (1988): 

Tb[/C] -1.22 X 10^2 (l + z)-A^ , (1) 

where S is the flux density in Jy, 9i and 02 are the sizes 
of the major and minor axes of the fitted Gaussian com- 
ponent in mas, and v is the observing frequency in GHz. 
These values are listed in Table 2. For the radio core. 
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Tb = 7.9 X 10^^ K, which is comparable to the inverse 
Compton Hmit (Kellermann & PauUny-Toth 1969). 

4. Discussion 

The overall radio structure of 3C 216 can be modelled 
with either helical flows or precessing beams in which 
jet components move ballistically, but a unique solution 
cannot be found (Fejes et al. 1992). Helical jets arc pro- 
posed to exist in some radio-loud AGN (e.g. Conway & 
Murphy 1993, and references therein) e.g. due to Kclvin- 
Helmholtz instabilities introduced by a driving mecha- 
nism such as precession of the central engine. Note that 
precessing jet models invoking ballistic ejection of plas- 
mons are less favoured at present as it would require two 
orders of magnitude higher density of the jet with respect 
to the ambient medium, which is not observed in AGN 
(Hardee et al. 1994). 

In the case of 3C216, Akujor et al. (1996) noted that 
the edge brightened structure and the flattened spectrum 
of the component near the jet bend may indicate inter- 
action with the interstellar medium. Such jet-ISM in- 
teraction is proposed to be a general phenomenon in the 
sample of CSS sources (including 3C 216) investigated by 
de Vries et al. (1999) with the Hubble Space Telescope, 
which gives evidence for the presence of emission-line gas 
aligned with the radio structures. The very high rota- 
tion measure observed in the jet bend (Venturi & Taylor 
1999) and the recently discovered H i absorption in the 
source (Pihlstrom et al. 1999) indicate that the jet is in- 
deed deflected by an interstellar cloud in the host galaxy. 
High resolution H i observations have been carried out in 
the UHF band with the EVN to identify the location of 
the absorbing gas (Y. Pihlstrom, priv. comm.). 



The position of component B from earlier VLBI mea- 
surements is shown in Fig. 3. These data are consistent 
with a constant proper motion of /U = 0.14±0.01 mas/yr 
between 1979 and 1999. But note that the higher res- 
olution observations (the two 8 GHz data points and 
our 5 GHz space VLBI measurement) show significantly 
larger separations of the component from the centre. The 
ground only data points obtained at 5 GHz are still in- 
dicative of jet deceleration or orientation changes, if we 
compare the trend in the component separation evolu- 
tion before and after 1988. It may be misleading to com- 
pare component separations measured with instruments 
of different resolution or measured at different frequen- 
cies, which may result in a shift of the reference point (the 
location of the observed brightness peak). This can hap- 
pen in a case e.g. when a jet component is blended with 
the core in a lower resolution observation. But the ref- 
erence point is also expected to shift with the observing 
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Fig. 3.. Motion of the supcrluminal component B from the earliest 
VLBI observations. The 1979-89 data arc taken from Barthel 
et al. (1988) and Venturi ct al. (1993), the 1991-98 data are 
taken from Venturi et al. (in prep.). Filled rectangles indicate 
observatons at 5 GHz, while asterisks at 8.4 GHz. 



frequency due to a change in synchrotron self-absorption 

opacity, if the core region can be described with com- 
pact jet models (Lobanov 1988 and referencies therein, 
see below). Such a shift was recently detected e.g. in 
the BL Lac objects 1803+784 (Perez-Torres et al. 2000) 
and 1823-F568 (Paragi et al. 2000). The latter was ob- 
served at four-frequencies (5-22 GHz) with the VLBA; 
these observations show a systematic increase of core-jet 
component separation with increasing frequency. 

The topic of constant component speed vs. jet decel- 
eration scenario for 3C216 will be investigated further 
by T. Venturi et al. (in preparation). For the present 
discussion we assumed the most simple case of constant 
component speed, and calculated the apparent transverse 
velocity with the following formula (e.g. Pearson & Zen- 
sus (1987)): 



Papp — M 



Ho{l + z) 



l + (l + 2goz)i/2 + z ' 
l + {l + 2qozy/^ + qoz. 



(2) 
(3.0± 



The observed proper motion corresponds to /?app 
0.2) apparent transverse velocity. 

The nature of component A is rather unclear. Its sep- 
aration from the core is 1.25 mas (see the super- resolved 
image in Fig. 2b) , compared to the earlier reported quasi- 
stationary position of 1.4-1.6 mas (Venturi et al. 1993). 
The beam orientation is not suitable to decide whether 
this feature is a separate component, or just an exten- 
sion to the core emission. As we could not model the 
source structure with a stable component at a distance 
of 1.25 mas in the model- fitting process, the latter seems 
more favourable. It is very probable that the "core" and 
"component A" together form the so-called ultracompact 
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Table 2. Fitted elliptical Gaussian model parameters of the source components. S: flux density, r: separation, Q: position 
angle, o: major axis, b/a: ratio of minor and major axes, $: position angle of the major axis, Tb: brightness temperature. 



Component 


S 


r 





a 


b/a 




Tb 




(mJy) 


(mas) 


n 


(mas) 




n 


(10^'' K) 


Core 


381.6 


0.00 




0.36 


0.31 


-32 


0.79 


A 


122.8 


0.87 


151.5 


1.37 


0.21 


-62 


0.03 


B 


45.1 


5.16 


152.3 


1.93 


0.14 


-31 


<0.01 



jet region of 3C216 (see Lobanov 1998). In this case we 
can apply the equipartition jet model of Blandford and 
Konigl (1979) in which the limiting brightness temper- 
ature is about 3 x 10" S'^^^ K. A comparison with the 
observed value results in a lower limit to the Doppler 
factor: 5 ~ 3.17. The Lorentz factor of the jet (7) and 
the viewing angle to the line of sight (9) can be given by 
(e.g. Daly et al. 1996): 

7 = (25)-i(/Jfpp + <52 + l), (3) 

and 

cos^=(-^) (7-r^). (4) 

Using the calculated values of /3app and the lower limit 
of S, and assuming h = 0.7 our results are 7 > 3.16 and 
< 18?4. Much higher Doppler factors for 3C216 have 
been determined by other methods {6 = 33 in Ghisellini 
et al. 1992 and ^ = 65 in Giiijosa & Daly 1996). These 
Doppler factors would result in very large Lorentz factors 
of 7 > 16 and much smaller angles to the line of sight 
(~ 1°), unless /3app represents a jet pattern speed instead 
of bulk motion. Note that the 5 estimations from the 
literature mentioned above made use of VLBI data taken 
during an outburst, when the flux density of the source 
core was much higher (cf. Barthel et al. 1988). 

5. Conclusions 

The blazar type core of the source can be interpreted 
in terms of compact jets. We estimate that the inner 
jet is oriented close to the line of sight with 6 < 18?4, 
and the radio emission is Doppler-boosted by a factor 
of ~ 3 at the epoch of observations. We conclude that 
the observed small projected size of 3C 216 is probably 
caused by both interaction of the radio beams with the 
interstellar medium and a projection effect. The nature 
of the core region should be investigated in more detail 
in future high resolution multi-frequency VLBI experi- 
ments, which would allow us to compare the observed 
spectral distribution with physical models. 
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